function [dx]=faii(~,x)
f1=x(2);
f2=-(x(1)^2-1)*x(2)-x(1);
dx=[f1;f2];
end